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Abstract. We propose a mechanism for fixing the velocity of relativistic 

■ soliton based on the breaking of the Lorentz symmetry of the sine-Gordon (SG) 
model. The proposal is first elaborated for a molecular chain model, as the 
simple pendulum limit of a double pendulums chain. It is then generalized to 
a full class of two-dimensional field theories of the sine-Gordon type. From a 
phcnomcnological point of view, the mechanism allows one to select the speed 
of a SG soliton just by tuning elastic couplings constants and kinematical pa- 
rameters. From a fundamental, field-theoretical point of view we show that the 

| characterizing features of relativistic SG solitons (existence of conserved topo- 

i logical charges and stability) may be still preserved even if the Lorentz symmetry 

is broken and a soliton of a given speed is selected. 



X 

H ' 1 Introduction 



In the past decades solitons have become a very useful concept in many areas of 
mathematical, theoretical and applied physics. Originally discovered as solitary 
waves solutions of non linear partial differential equations [TJ [51 [3] , solitons have 
become a generic paradigm of non linear physical phenomena [H [5] . 

Of particular relevance, both from the theoretical and phcnomcnological 
point of view, are relativistic solitons, in particular topological two-dimensional 
(2D) sine-Gordon (SG) solitons (kinks) [6], on which this paper will be mainly 
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focused. On the phenomenological side, SG solitons have wide range of applica- 
bility: non linear molecular [Jj and DNA dynamics [H[9l[T0l[Tl], Josephson effect 
[T2"] , ferromagnetic waves [T31 [H] , non linear optics [T51 HE] , superconductivity 
[T7] and many others [TS] . 

From a theoretical point of view SG solitons exhibit remarkable analogies 
with relativistic particles. They share with the latter the feature of being lo- 
calized, non dispersive objects realizing the Lorentz symmetry. SG solitons are 
characterized by short range forces and they can interact without loosing their 
identities [19] . 

Moreover, topological SG solitons are stable and their existence is related to 
the discrete internal symmetries (and the existence of degenerate vacua) of the 
two-dimensional (2D) field theory of which they are classical solutions. Topology 
offers an elegant way of classifying solitons in terms of the homotopy group of 
the mappings between the vacua and the asymptotical field configurations [6]. 

The soliton-particle analogy has been very fruitful both for particle physics 
and quantum field theory. Starting with the pioneering work of Finkelstein and 
Misner [5D] , Rubinstein \§\ and Skyrme [3T] , the analogy has found applications 
in various contexts such as S-Matrix formalism [22], current algebra effective 
Lagrangians [23] . supersymmetric quantum mechanics |24j . black hole physics 
[231 12"5] just to mention a few of them. 

Solitons have also became central in the study of integrable PDEs [271 HE] 
and in many geometrical matters (in particular connected with Field Theory), 
see e.g. [29] . 

A general feature of relativistic solitons and in particular of SG solitons, is 
the fact that the soliton speed v is a free parameter, which can be fixed by 
choosing initial conditions and is bounded from above by a limiting value Co- 
This is obvious in view of the Lorentz symmetry of the model: only the limiting 
upper bound (corresponding to the speed of light for a relativistic particle) of 
the soliton speed is fixed by the model, whereas v can be changed applying a 
boost. 

This degeneracy in v is a necessary property of a relativistic field theory 
such as SG theory and appears rather natural in view of the soliton-particle 
analogy. However, in phenomenological situations where the experiments give a 
well-defined value for the propagation speed of the soliton it represents a loss of 
predictive power of the model. By modelling some non- linear dynamical process 
using a SG-like field theory very often it is crucial not only that solitary, non 
dispersive waves do exist but also that their velocity is at least of the same order 
of magnitude of that observed in experiments. 

Thus, a selection mechanism for the soliton velocity should be very welcome 
from the phenomenological point of view. In this context an instructive example 
is represented by the use of solitons in the description of DNA non linear torsion 
dynamics. It is believed that kink solitons of the SG type play a crucial role 
in the DNA transcription process [HI IS] • In the past decades several models of 
DNA torsional dynamics have been proposed that allow for the existence of SG 
kink solitons ( see e.g [HUH] and references therein). However, in most of them 
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the soliton velocity is not fixed by the mode£j], whereas the DNA transcription 
process, which the soliton propagation is assumed to model, occurs with a speed 
of a well-defined order of magnitude. 

Although desirable from a phenomenological point of view, it is not very 
easy to implement a general selection mechanism for the soliton speed. The 
high sensibility of non linear systems to the initial conditions, practically rules 
out a selection mechanism based entirely on tuning initial conditions of the 
system. 

On the other hand, one expects a speed selection mechanism to have a strong 
impact on the existence and on fundamental properties of SG solitons. In fact 
any velocity selection breaks the Lorentz symmetry of the model. It is not a 
priori evident that after this breaking the theory will still allow for localized 
solutions with the features of the SG solitons such as stability and topological 
classification. 

In this paper we will present and discuss in detail a velocity selection mech- 
anism for SG topological solitons. Our starting point is a result we found in 
previous investigation of a composite model for DNA torsional dynamics, which 
is essentially a double chain of coupled pendulums [301 OS] ■ We found that in 
the limiting case where one of the two coordinates describing the angular dis- 
placement of the pendulums is constrained to be identically zero, the continuum 
limit of the model has SG solitonic solutions with speed fixed in terms of the 
elastic and inertial parameters of the chain. 

The possibility of selecting the speed of SG solitons by adding to the La- 
grangian terms breaking the Lorentz symmetry has been first recognized in Ref. 
[32] , A speed selection mechanism has been also found in the case of solitonic 
solutions of the non-linear Schrodinger equation [33j . In particular, in Ref. [32] 
the authors consider a SG model modified by the addition of further kinetic and 
potential terms. In this paper we consider instead a two fields generalization of 
the usual SG model. Moreover the breaking of the Lorentz symmetry, which is 
responsible for the fixing of the soliton speed is rather subtle. It is generated 
by the presence of two (instead of one) upper limit speeds for the propagation 
of waves. 

The plan of the paper is as follows. We will first review the results of Ref. 
[30] and reformulate them as the simple pendulum limit of a double pendulums 
molecular chain model (Section [2]). We will show that the selection mechanism 
is related to the existence of a conditionally conserved quantity, a conservation 
law that holds only when the system is constrained in a limited region of its 
phase space, i.e. to a region of the reduced one-dimensional (ID) equations of 
motion describing propagating waves (Section [3]). We will then pass to consider 
the most general model for which the selection speed mechanism does work. We 
will show that the general model can be formulated as a two fields generalization 
of 2D SG field theory (SectionUJ). We will consider a simplified, minimal version 
of the model, which has the nice feature of being explicitly solvable in its general 

1 The only exception is the model of Ref. [30] , where the selection mechanism discussed in 
this paper was first observed but not identified in its generality. 
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form (Section [5]). The breaking of the Lorentz symmetry is discussed and it is 
shown that it is related to the presence of two (instead of one) upper bound 
limiting speeds (Section 0. The discrete symmetries of the model are discussed 
and used for the topological classification of the kinks (Section [7]). The stability 
of our solitonic solutions are also investigate. In particular we show that the 
Lorentz symmetry breaking term acts as damping term for linear perturbation 
near the soliton (Section [8]). Finally we state our conclusions (Section [9]). 

2 Soliton speed selection in molecular chain 
models 

Let us start from the double pendulums chain model of Ref. [30] modeling DNA 
torsional dynamics. The two DNA helices are modeled by two chains of double 
pendulums. At each, equally spaced with distance 5, site of the two chains 
there is a first pendulum, which is a disk radius R and momentum of inertia I. 
Attached to the border of the disk there is second pendulum of length r and mass 
m. The rotation angles of the two pendulums with respect to the equilibrium 
positions are denoted respectively by 9f, ip", where the index a = 1,2 refers 
to elements in the two chains and i £ Z identifies the site in the chain and 
Of,<pf £ S . The pendulums interact with nearest-neighbor interactions of three 
types: stacking (coupling between successive pendulums on the same chain), 
pairing (coupling between pendulums on the same site but opposite chain) and 
torsional (coupling between successive disks on the same chain), characterized 
by coupling (elastic) constants given respectively by K s , K p , K t . Moreover, the 
second pendulum is not free to swing through a full circle but it is instead 
constrained to stay in a range tpf < ipa < 2n. This constraint can be modelled 
by adding a confining potential V c (<p), which has the effect of limiting de facto 
the excursion of the <p angles; see Ref. [30J |3T] for details. 

As here we are interested not in the details of the model but in its general 
behavior, we will make, with respect to Ref. [301131] , a number of simplifying as- 
sumptions. First we will consider from the beginning only symmetric solutions, 
0} = 9f := $i, ip\ — <pf := ifi. This position reduces the degrees of freedom 
from four to two per site. Second, to reduce the number of kinematical param- 
eters we will set R = r. Under these simplifying assumptions, in the continuum 
(S — > 0) limit the dynamics of the model is described by a 2D field theory with 
Lagrangian density 

c = ^{ie^-Lu t e 2 x + r 2 [me 2 t -uj s e 2 x + 2cos(ip)(me t (ipt + e t ) + 

- uj s 9 x ((p x + 6 X )) + m(ip t + 9 t ) 2 - Lo s {tPx + X ) 2 ] } + 

+ Ar 2 K p ^cos<9 + cos(<p + (9) - ^cos(ip) - ~)^ -V c (<p), (1) 

x is the spatial coordinate along the chain, the fields 9(t, x) and tp(t, x) represent 
the angular coordinate 9i and <pi in the continuum limit, the indices x and 
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t denote partial derivation and uj t = K t S 2 , u> s = K s 6 2 . We will use in the 
Lagrangian ([1]) a rather generic confining potential V c . We will only assume 
that V c (ip) is an even function of tp, which has a minimum at tp = 0. Moreover, 
we can take without loss of generality 14(0) = 

At this point the interpretation of the Lagrangian |T|) as modelling the DNA 
torsional dynamics is not anymore compulsory, nor so relevant. We can as well 
consider the Lagrangian as describing the continuum limit of the torsional dy- 
namics of a single molecular chain made of two pendulums. The pairing interac- 
tion for the DNA double chain becomes an external potential for the pendulums: 
V = -~4r 2 Kp [cos 8 + cos(</j + 8) — i costp — 3/2), whereas the pairing and tor- 
sional interaction generate in the Lagrangian ([1]) the gradient (x-derivative) 
terms. 

Let us consider travelling waves solutions of the field theory (JTJ) , i.e solutions 
depending only on z = x + vt with fixed speed v: 8 = 8(z), tp = tp(z). With this 
position the field equations stemming from {!]) become 

flip" + /Lt(l + cos tp) 8" = 

= -AK p sin (tp + 6) — fi sin (tp) (8'f + 2K p sm{tp) - ; 
H(l + cos tp) tp" + [(J / r 2 ) + 2 t i(l + cos tp)} 8" = 1 ' 

= -AK p (sm8 + sm(tp + 8)) + fi sm(tp)[(tp') 2 + 2tp'9'] , 

where the prime denotes derivatives with respect to z and following the notation 
of Ref. [315] we define 

fi := {mv 2 —u}„), J := [Iv 2 - u t ) . (3) 

In general the system ([2]) cannot be solved analytically in closed form. One 
has to resort to numerical calculations in order to show that the system admits 
solitonic solutions [30]. However, in Ref. [30] it was shown that one can find a 
solitonic solution with fixed speed, just by freezing the angle tp, i.e by setting 
tp = 0. We note that if we force tp(z) = 0, we are actually considering a 
chain of simple pendulums, i.e. a sine-Gordon equation. This constraint can be 
accommodated in our setting in a dynamical way, by acting on the confining 
potential V c : this should be made stronger and stronger and the maximum angle 
tpo smaller and smaller. In the limit tpo — > + and (d 2 V c /dtp 2 )(0) — ► +oo, we 
expect to recover the solitons of the sine-Gordon equation. 

Setting tp = and using (dV c /dtp)(0) = 0, the system ([2]) is equivalent to 
the SG field equation 

fiB" = -2K P sin 9, (4) 

with the compatibility condition between the two equations of the system ([2]) 
given by 

J = 0. (5) 

In principle one could consider a generic constraint tp = q8, where q is a pro- 
portionality constant. However, it is not difficult to realize that for q =/= such 
a constraint, when used in Eqs. |[2]), will lead to two different determinations of 
8 making the two resulting equations incompatible. 
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Eq. Q has to be integrated with the usual boundary conditions for kinks 
6*(±co) — 2u±tt, 6*'(±oo) = with n± 6 Z. Here n = n + — n_ is the kink 
winding number. For n = 1 Eq. (|U) can be easily integrated to give the kink 



6> = 4arctan[e' Cz ], 



9? = 0, K- 




(G) 



whereas the compatibility condition ([5]) fixes the speed of propagation of the 
soliton to the speed Ct of the transverse sound waves supported by the elastic 
torsional forces acting on the disk 



Moreover, because the soliton exist only for /i < 0, the soliton speed is bounded 
from above by the speed c s of transverse sound waves supported by the elastic 
stacking forces acting on the pendulum, i.e. 



In view of Eq. (JT]) , this implies the following constraint on the stacking, torsional 
coupling constants and kinematical parameters of the system: 



The condition (p — fixes the speeds of the soliton to the value (|7|) also 
for kinks with winding numbers n > 1. In fact, in Sect. [7] we will show that 
the fixed speed solitonic solutions of Eq. (U|) allow for the same topological 
classification of the usual SG solitons. 

The selection mechanism for the soliton speed discussed above represents a 
very nice and simple way to produce SG solitons with a given speed in double 
pendulums molecular chains. To select the soliton speed one just needs to tune 
the torsional and stacking coupling constants and the kinematical parameters 
of the chain such that Eqs. (O and ([9]) are satisfied. Acting on the confining 
potential V c , making it stronger and stronger, one obtains the single pendulum 
limit of the double pendulums chain. The angle is frozen to if = and a SG 
soliton with a speed equal to that of the transverse sound waves supported by the 
torsional forces acting on the disk is selected. Notice that the above described 
mechanism can be obviously used to devise and realize non-linear media where 
solitons propagate at a given fixed speed. 

Although very simple the mechanism we are proposing raises a number of 
questions. How is our mechanism related to fundamental features of the SG 
system such as integrability, Lorentz symmetry and topological classification? 
Can our mechanism be generalized to a broad class of models or is it just an 
artifact of the peculiar couplings present in the Lagrangian ([1])? Is the SG 
soliton we are selecting stable? We will try to answer these questions in the 
following sections. 



V = Ct = \fuhfl. 



(7) 



v < c s = y/uJs/rn. 



(8) 
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3 Integrability and conditionally conserved 
quantities 

The simple (one single degree of freedom) SG system ^ is integrable, in the 
sense that the initial value problem can always be solved analytically in closed 
form. Moreover, this integrability is related to the existence of a conserved quan- 
tity (a first integral of the equations of motion) , which can be easily identified 
as the energy of the system. 

On the other hand, the (two degrees of freedom) system in general allows 
for only one conserved quantity - the energy - and it is not integrable in the 
sense used above. Nonetheless, it can be solved analytically in closed form 
and reduces to the SG system (0| if we limit ourselves to consider only tp = 
solutions. One is therefore led to search for a conditionally conserved quantity 
(a discussion of this weak form of integrability can be found in Rcf. 34 ) . That 
is, a quantity that in general is not conserved by the dynamical evolution of the 
system, but which becomes conserved when the dynamics is restricted to the 
subspace tp = 0. 

To this aim, let us first note that the equations of motion ||3J| describe a 
dynamical system with two degrees of freedom, with the coordinate z playing 
the role of the time. The equation of motion can be derived from the Lagrangian 
L = T — V, 

(6') 2 + 2cos(^) ((6 1 ) 2 + 6'tp') + (9' + p'f 
- 2K p r 2 [costp - 2cos(6> + tp) - 2cos6»] - V c {tp). (10) 
With a little algebra we obtain from the system @ the equation, 

^P(0^^)=f(^,(fl') 2 ,^') 2 )- (ii) 

where J- is given by 

dV 

T = 2r 2 [2K p (sin(0 + tp) - sin 6 - sin tp) + fismp((9') 2 + tp'6')] +2-rp. (12) 

and P is linear in the momenta Pg = 8L/S6', P v — SL/Stp': 

P = P e -2P v = J9' + s(ip)fiip', s = r 2 {costp-l). (13) 

Notice that the coefficient s, its derivative with respect to tp and T vanish when 
evaluated on tp = 0. The physical interpretation of Eq. (|11[) is very simple: this 
is just Newton's second law with P playing the role of the impulse and T that 
of the force. 

The equations of motion ^ can be now rewritten in the equivalent form 

fitp" + /i(l + cos tp)9" = -TftTpSin (tp + 0) - fisin (tp) (9 1 ) 2 + 2K p sm(tp) 

P' = T. (14) 



l = - n J(e'? + ^ 
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In general the momentum P is not conserved, but becomes conserved when 
T = 0. This can be achieved if we restrict ourself to consider only the solutions of 
Eqs. that satisfy tp = 0. Taking into account that F(0) = s(0) = s'(0) = 
Eq. (fTTj) gives a conservation law for 

P(ip = 0):=H= JO'. (15) 

Thus, H is the conditionally conserved quantity we have been searching for. 

It can be interpreted has an "effective" angular momentum of the disk. In 
fact J has two contributions, see Eqs. ©. The first comes from the geometrical 
momentum of inertia of the disk, the second has opposite sign and comes from 
K t , the elastic forces acting on the disk. 

For ip — 0, Eqs. flip reduce to 

dU 

fi6" = -2K P sin 9, — = 0. (16) 

For J ^ the previous equations do not admit solitonic solutions (actually the 
only solution is the trivial one 8 = 0). Solitonic solutions exist only for J = 0, 
when the conservation law is trivialized. In this case Eqs. (|16p coincide with 
Eq. ([4]), ([5]), which give the fixed speed soliton solution ([6]). 

Thus the physical mechanism behind the fixing of the soliton speed is the 
following. When restricted to the subspace of solutions <p = the system 
has a conservation law, which is not compatible with the existence of solitonic 
excitations. The requirement of existence of solitons implies a trivialization of 
the conservation equation, which in turn fixes the soliton speed. From this point 
of view the condition J = can be seen as a sort of resonance condition: both 
the force and the "effective" mass in equation (fTTj) must be zero. 

4 The general model as a 2D SG-like field theory 

The mechanism for selecting the soliton speed described in the previous sections 
for the molecular chain model UJ is rather generic. As we have shown in the 
previous section, the speed selecting mechanism is related to the existence of a 
conditionally conserved quantity and it is rather independent from the specific 
form of the interactions characterizing the model. Therefore, it will be present 
if the system satisfies some general conditions. 

4.1 The general speed selection mechanism in molecular 
chains 

Below is a list of the main ingredients that are needed for the mechanism to 
work in the case of a molecular chain. We need a chain (array) of identical 
mechanical systems with the following features: 

1. The system must have at least two degrees of freedom (X, Y) characterized 
by two masses (or momenta of inertia) m, M with m ^ M; 
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2. There are at least two types of interactions: a) An elastic force (coupling- 
constant K t ) originated by the interaction between neighboring sites on 
the chain; b) A non linear external force (coupling constant K p ) acting on 
the single site; 

3. There is a confining potential V c that limits the range of variation of one 
degree of freedom (e.g Y) and allows to freeze Y. That is by making V c 
stronger and stronger we perform the dynamical reduction from two to 
one single degree of freedom X: 

4. Freezing the degree of freedom Y we have both a conservation law for the 
momenta conjugate to X and solitonic solutions for X . 

If the conditions above are satisfied, then the mechanism described in the 
previous sections will work. In particular, considering the continuum limit of 
the chain and travelling wave solutions X(z) — X(x + vt), Y(z) = Y{x + vt) we 
get a system of two differential equations which describe a mechanical system 
with two degrees of freedom with effective masses d 

m cff = v 2 m - K t S 2 , M cff = v 2 M - K t 8 2 . (17) 

As the conservation of the momentum conjugate to X (for Y frozen), Px — 
M c f[dX/dz is not compatible with the existence of solitons, this will require 
M e & = 0, which will fix the soliton speed in terms of transverse phonon speed 
of the chain, 




Moreover, m e ff is the analogous of the parameter \i in Eqs. ([3]) and {!]), hence 
the soliton will exist only for m e g < 0, implying v < 5^jK t /m. Consistency of 
these equations requires M > m. 



4.2 Speed selection in field theory 

The explicit construction of models satisfying the previous conditions becomes 
much simpler if we consider them directly at the field theoretical, Lagrangian 
level, as generalization of the SG system, rather then as originated from a molec- 
ular chain. The main advantage of this approach is that the speed selection 
mechanism becomes a general, field theoretical, effect not necessarily originated 
from the continuum limit of some discrete model. In the rest of this paper we 
will follow this approach and we will make no more reference to molecular chain 
models. 

To write down the Lagrangian for a general class of field theories that exhibits 
the speed selection effect, we first redefine the field ip in the Lagrangian density 

2 Depending on the particular model one is considering, the equations below may also 
assume a more general form m e g = v 2 m — f(ri)Kti 2 (and similarly for M e fj), where /(r;) is 
a function of the geometrical parameters r% characterizing the model. 
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|T]). Defining cp = $ — 9 the Lagrangian |T]) becomes 
C = - {l9 2 - cj t 6l + r 2 [m6 2 - oj s 9 2 x + 2cos($ - 9) (mO t $t - 



1 3 

-m<& 2 - lu s $1 + 8K p ( cos 6 + cos $ - - cos($ - 6) — - 



-V c {*-0). (19) 



In this new parametrization, the fixed speed soliton solution is given by $ = 
9 = 6» , with 9 Q given by Eq. ©. 

Solitons are localized excitations of finite energy. The energy of the solution 
([6]) can be calculated from the Hamiltonian density. From the Lagrangian (| 19|) 
one can easily obtain the Hamiltonian density 

H = - {I6 2 t + oj t 9 2 x + r 2 [m£ 2 + oj s 6 2 x + 2cos($ - 9) {m6 t $ t + uj s 6 x $ x ) + 



+ m$ 2 + u>s$ x - 8K P ( cos 9 + cos $ - - cos(<E> 



1 N 3 

-cos(^-e)-- 



+ V c ($-9). (20) 



The constant field configuration of minimal energy (vacua) can be obtained by 
minimizing the potential term in Eq. (|20[) . They are given by 9 = $ = 2m:, n 6 
Z and have zero energy. The energy of the soliton ([6]) , connecting vacua with 
n+ = ri- + 1, is 

H= r Hdx= 16 r 2 K p ^ + ^ , (2 l) 
J-oo u s - {m/I)uJt 

As expected the energy of the soliton (measured with respect to that of the 
vacuum) is finite and positive in the range of existence of the soliton. 

It is useful, in particular in view of the discussion about the Lorentz symme- 
try of the theory, to write the Lagrangian using a Minkowski spacetime nota- 
tion. Introducing a spacetime metric of signature (1,-1), the two fundamental 
speeds c t , c s (we assume c t < c s ) given by Eqs © and respectively ([7]), and 
the derivative operators 

the Lagrangian (fTTJ|) takes the form 

2 

£ = ^-d v 6d v 6 + {d u 9d u 9 + 2cos($ - 9)d u 9d v <S> + d„$d v $) 



+Ar 2 K p ( cos 9 + cos $ - t CO s($ - 6) - ~ ) - F c ($ - 9). (23) 



Generalizing the Lagrangian (f2"3"|) one obtains a class of models admitting soli- 
tonic solutions with fixed speed. 

The models have the form of a sine-Gordon-like coupled system of two scalar 
fields and will be functionally parametrized by two coupling functions F, G. In 



10 



the Minkowski spacetime notation the Lagrangian is given by 

C g = ~d v 6d v 6 + h -d v Qd v e + F($ - d)d v ed v <S> + 

+ ~<9„$d"$ + is:(cos0 + cos$-2) + G($-0), (24) 

where a, 6, K are positive constant parameters and the coupling functions _F($ — 
9) and G(«J> — 6*) are arbitrary even functions of <& — 6 and G satisfies G(0) = 0. 
Considering only travelling wave solutions <&(z), 0{z) (with z — x + vt) we obtain 
the ID reduced Lagrangian which is the generalization of the one given in (fT0|) : 

L = -J{d'f + | [{9' f + 2F6'®' + ($') 2 ] + K [cos $ + cos(0) - 2] + G, (25) 
where 

J = a(v 2 /c 2 t - 1) , £ = fe(« 2 /^ - 1) • 

Taking into account that because of the parity of the coupling functions F and 
G we have F'(0) = G'(0) = 0, it is not difficult to show that for $ = 6 the 
equation of motion stemming from the action (|2"5|) admits the SG soliton ^ 
travelling at speed v — c t < c s . 

A further generalization, which we will not investigate in detail in this paper, 
can be obtained considering generically a self-interacting scalar field 8 in 2D 
Minkowski spacetime, 

c = -d v ea u e - v(e), (26) 

whose potential V is such that field equations for O admit relativistic soliton 
solutions (apart from the SG system, the model contains as a particular case 
e.g the ip* interaction model) . A general Lagrangian describing two interacting 
scalar fields 0, S for which the speed selection mechanism for the soliton applies 
is given by 

-V(Q) - V(E) - G(Q - E) , " 1 ' 

where as above F(Q — E) and G(0 — E) are arbitrary even functions of (9 — E). 

5 Minimal model 

The Lagrangian (|24|) corresponds to a system of two coupled sine-Gordon mod- 
els. It is of interest to see if the properties concerning the selection of the 
soliton speed hold also when the interaction terms are switched off, i.e when 
the coupling functions F, G are identically zero. Moreover in the completely 
non-interacting case, the system will be integrable and we will have complete 
control on the general solutions. 
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For F = G = the Lagrangian (jiM)) becomes 



C m = ^b v 9d v 9 + h -d v dd v e + -d v §d v ® + K (cos6 + cos $ - 2) . (28) 



The field equations describe a decoupled system for 9 and <&, 



-d v d v e + d v d v o 



K 



sine, ft,a"* 



A" 



• sin 



(29) 



The theory contains two upper bound speeds, c s ,c m for the propagation of 
respectively $- and #-waves. c m is given by 



CtC s 



act 



bcT 



(30) 



Notice that c s > c t implies necessarily also c s > c m . 

Because the two scalars are decoupled we can look for traveling wave solu- 
tions of different speeds v, v for the two fields: $(z) = &(x+vt), 9(z) = 9(x+vt). 
In this way we we end up with a decoupled system of differential equations, 



« \ n K 

-j - 1 *" = -— sin$, 

ct o 



- 1 



- 1 



K 

= — — sin 
b 



9, (31) 



where the prime and the dot denote derivation with respect to z and z respec- 
tively. The general solution of the previous system is given by two sine-Gordon 
solitons of type ^ with speeds v < c s and v < c rn < c s , respectively for the 
fields $ and 9. Only imposing 9 — $ we have velocity selection. In fact 9 — $ 
implies v = v = c t Notice that the converse is not true, v — c t does not neces- 
sarily imply 9 — $ or v = v. For v = c t and 9 ^ $ we will have solitons in 9 
and $ that differ in speed or, if the speed is the same, in the coordinate of the 
center of the soliton ( we will have two identical solitons with non-zero relative 
phase). 

6 Lorentz Symmetry 

It is well known that usual SG field theory, described by the 2D action 
'(ft) 2 1 



S 



2cl 2 



(ip x ) +w (cosy>-l) 



dxdt , 



(32) 



is invariant under 2D Lorentz transformations (boosts) 



x=j(x' + vt'), t = j(t' + —x'), 7 



7 1 

°0 



(33) 



where cq plays the role of the speed of light in vacuum. The Lorentz symmetry 
of the action (|32|) is essential for establishing the soliton/relativistic particle 
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analogy and allows one to use for SG solitons the concepts of relativistic kine- 
matics. For instance, Lorentz symmetry allows us to generate, using a boost 
(|3"3"|) , a SG soliton propagating with speed v 

<p(x + vt)=4 arctan [^7(^+^)1 ) (34) 

from the static solution ipo(x) — 4arctan[e" o2: ]. Another consequence of the 
Lorentz symmetry of the theory is the existence of a maximum propagation 
speed, v m ax = Co for the solitons. To make the Lorentz symmetry more evident 
we might write the action (|3^|) using the Minkowskian notation of Eqs. l|2"2"]) . 

On the other hand, the two-fields generalized SG theory we are dis- 
cussing in this paper is not Lorentz invariant. This is due the fact that in the 
Lagrangian (j2"4")l appear not a single upper bound speed Co but two of them 
c s ,c m . It follows that we have the choice of defining the boosts (|33[) either in 
terms of c s or c m . In the first case the Lorentz symmetry is broken by the kinetic 
term proportional to a, in the second case by the kinetic terms proportional to 
b in the action . Because c m < c s it is natural to consider boosts of the first 
kind, so that it is the o-term that breaks the Lorentz symmetry. 

From the physical point of view, the breaking of the Lorentz symmetry can 
be easily understood if one realizes that our system can be thought as a medium 
with two different speeds for the propagation of sound waves. 

It is obvious that the speed selection effect described in previous sections is 
strongly related to breaking of the Lorentz symmetry. In fact if in Eq. (|24[) 
we set a = we have no selection for the speed of the soliton and Lorentz 
symmetry is not broken. Let us therefore investigate in detail the relationship 
between breaking of Lorentz symmetry and soliton speed selection mechanism, 
starting from the minimal model 



6.1 The minimal case 

In the minimal case, the fields, $ and 9 do not interact one with the other. This 
means that the field theory has two completely decoupled sectors. The field 
equations ([29]) can be written as follows, 

1 K 1 K 

-irO tt -O xx = —r sin9 > -»®tt-®xx = —r sin$ > ( 35 ) 

cfn a + b cj b 

where c m is as in Eq. f|30[) . Notice that we have three fundamental speeds in 
the theory, viz. c m ,Ct,c s with c m < c t < c s . Here c m and c s are the maximal 
speeds of propagation respectively for 9- and <I>-waves. The physical meaning of 
Ct will be clear later. 

Because of the decoupling, the two sectors of the field theory are invariant 
under two different groups of Lorentz transformations. Using a notation in 
which the dependence of a boost from the upper bound speed is evident, we 
will denote with r(c ) the boost in (|3"3"|) . Whereas the field equation (|3"5|) for 
9 is invariant under boosts r(c m ), that for <!> is invariant under r(c s ). The 
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previous symmetry properties have a simple physical consequence: the existence 
of solitonic excitation for the fields 9 and $ which propagate at two independent 
speeds v, v (see section[5]). The $-soliton, $>(x+vt) can be obtained boosting the 
static soliton <&(&) with r(c s ), whereas the fl-soliton, 9(x + vt) can be obtained 
boosting the static soliton 9(x) with T(c m ). The solutions can be written in an 
Lorentz invariant form introducing the covariant vectors 

k^ = (k,u), y i = (k,Q), x^ = {c s t,x), x^ = {c m t,x). (36) 

Propagating solutions are function of the Lorentz scalars, 

$ = $(k^,k^) , e = 6(k^, k^) . 

The condition 9 = <!> couples the two otherwise decoupled sectors and breaks 
Lorentz invariance. In fact 9 — $ requires v = v, whereas Lorentz invariancc 
would require r(c s ) = P(c m ), which is manifestly impossible for c s ^ c m . As 
shown in Section [5j the condition 9 = $ fixes the speed of the soliton so that 
we have v = v = c t . The soliton propagating at this speed cannot be obtained 
by boosting a static solution. This is also evident considering that Eqs. (|3"Tj) do 
not allow for static solutions with <!> = 9. 

From a physical point of view the mechanism described above can be seen 
as a speed selection effect generated by breaking of the Lorentz symmetry. We 
have a media with two different maximal propagation speeds, c s , c m that can 
be identified as the propagation speed of $- and #-sound waves. Furthermore, 
we have a constraint, $ = 9 which breaks the Lorentz symmetry and singles out 
a solitonic perturbation which can propagate only at the fixed speed c*. 

6.2 Interacting models 

Until now we have considered only the minimal, non interacting model (|28[) . 
Passing to consider the interacting model (|24[) there is only one main difference: 
as now the fields 9 are coupled to each other, only one single boost T(c s ) can 
be defined, and the Lorentz symmetry is broken already at level of the action 
Being mutually coupled, travelling wave solutions for 0,$ now have to 
propagate with the same velocity v. Using the covariant notation of Eq. (|36[) 
they take the form 

$ = k^, uj 2 c 2 s /cI - k 2 ) , 9 = 9(k^,k^, uj 2 c 2 s /c 2 m - k 2 ). (37) 

Owing to the dependence on the non-Lorentz invariant term (u} 2 c 2 /c 2 n — k 2 ) 
the solutions are not Lorentz- invariant. Although the boosts r(c s ) are not 
a symmetry of the Lagrangian ([M]) they can be used as solution-generating 
transformations. A solitonic solution travelling with speed v can be obtained 
boosting a static solution. Similarly to what happens in the minimal case, the 
condition 9 = $ breaks this solution-generating symmetry. The 9 — $ solution 
cannot be generated by boosting a static solution satisfying the constraint, so 
that the velocity v = Ct for the soliton is selected. 
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7 Discrete symmetries and topological classifi- 
cation 



The internal (discrete) symmetries of our model are very important both for 
understanding the peculiarities of the fixed-speed 9 = $ solution and for the 
topological classification of the solitons. 

We are dealing with a 2D field theory with broken Lorentz symmetry, for 
which the internal symmetry group of the Lagrangian (|24[) is not necessarily the 
same as that of the reduced action (|25[) describing travelling wave solutions. We 
are only interested in the symmetries relevant for the solitonic solutions, namely 
in the symmetries of the reduced Lagrangian (|25[) . 

For J / and for generic coupling functions F, G the internal symmetry 
group of the model - denoted as Q in the following - and the residual symmetry 
of the vacua - denoted as Ti below - are the same as those of the simple SG 
system (|32|) . We have Q = Z x Z2, where Z is realized as discrete translations 
$ — * $ + 27m, 9 — ► 9 + 2nn, whereas Z 2 is the inversion 9 — » —6, $ — * The 
vacua $ = 9 = 2mir are invariant under an inversion followed by a translation 
of 47rm. The residual symmetry group of the vacua is therefore Ti = Z2. 

For J = the group Q acquires an additional Z2 factor: Q — Z x Z2 x Z2. 
This third factor comes from the invariance of the Lagrangian ([25]) under the 
transformation 

e < — ► $ . (38) 

The group Ti acquires an additional factor as well: Ti = Z2 x Z2. Notice 
that the speed-selecting constraint leaves invariant both Q and Ti. In fact the 
submanifold $ = 9 is selected by the constraint. 

The solitonic solutions of the model can be therefore classified in two classes, 
which are mapped one into the other by the transformation (|38|) . The solitonic 
solutions with 9 — $ are self-dual, i.e a fixed point of the transformation (|3"5|) . 

The previous discussion holds true for generic coupling functions F, G. For 
particular choices of F and G we can have bigger Q and Ti groups. This happens 
for instance when F and G are periodic functions. In this case Q is the product 
of two Z groups times some Z2 factors. However, the 9 = 4> condition breaks 
invariance under one of the two Z groups and we are left with the same symmetry 
groups Q and Ti as above. 

The discrete symmetries of SG-like field theories are also of fundamental 
importance for the topological classification of the solitonic excitations of the 
theory. The admissible number of solitons of the SG model is determined by the 
number of ways in which the points x = ±00 can be mapped into the different 
vacua of the model. 

We have therefore a one-to-one correspondence between solitons and the ele- 
ments of the homotopy group TTo(Q/Ti). It is also well-known that field theories 
that admit topological solitonic solutions have conserved currents correspond- 
ing to conserved topological charges, which can be identified with the winding 
numbers of the soliton. 
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Using the previous results one can easily calculate the homotopy group 7To 
in the case of generic coupling functions F, G. For J ^ this is given by 

MG/n) = 7T (z x z 2 )/z 2 = z . 

The same result holds true when J ^ 0: 

MG/H) = ir (Z xZ 2 x Z 2 )/(Z 2 x Z 2 ) = Z . 

As expected the solitonic solution of the model can be labeled by an integer n, 
its topological charge. 

For periodic coupling functions the homotopy group ttq of the model becomes 
Z x Z and we have two topological charges (n, m) associated to the soliton. 
However, the condition 9 = $ reduces the symmetry group and we are left 
again with ttq = Z and a single topological charge n. 

8 Stability 

The stability of the solutions we have considered in the previous section is a cru- 
cial requirement for our speed selection mechanism. First, being the dynamics 
non linear we expect high sensitivity of the solutions to the initial conditions. 
If the 8 — $ solution is not stable, small perturbations in the initial conditions 
may disrupt the mechanism. Second, even though a soliton with the given speed 
is generated, it may be unstable under perturbation and decay in a finite time. 

In the previous section we have seen that the fixed-speed solitonic SG soli- 
ton (|6|) , although resulting from breaking of the Lorentz symmetry, share many 
feature with the usual SG solitons: they are localized, finite energy solutions 
of the same ordinary differential equation, they allow for the same topological 
classification and associated topological charges. These features strongly sug- 
gest that the fixed-speed SG solitons we are discussing in this paper are stable 
solutions, in the same way as relativistic SG solitons are. 

The analysis of the stability of the solutions of the model (|24p is rather 
complicate. In this paper we will restrict our investigations to the classical 
stability of the speed-selecting solution 6 = $ under linear perturbations. To 
begin with, we will briefly summarize well-known results about the classical 
stability of the Solitons of the simple SG model (|32l) . 

To discuss the stability one first uses Lorentz invariance to chose the reference 
frame in which the soliton (|34[) is static and then considers small perturbations 
around the static solution, 

<p(x,t) = <p (x) + ip(x,t), ip « 1. (39) 

Inserting this position in the field equations, considering only the linear approx- 
imation and separating the variables in the perturbation, i.e setting ip(x, t) = 
\& '(x)e luJt , one obtains an equation for \t, which has the form of a time- inde- 
pendent Schrodinger equation 

/ m\ 2 d 2 
Tq,= \7) *' T ^ + UJ o( 1 - 2s ^h 2 x). (40) 
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Classical stability requires the eigenvalues of the operator T to be real and non- 
negative. This is guaranteed if the operator T is self-adjoint and non-negative. 
These two properties of the operator T given in Eq. (|4H|) can be shown in 
different ways. 

The most elegant way is to use supersymmetric (SUSY) factorization [2H 
[551 131)] . One can find a superpotential W(x) such that T = A + A with A = 
(d/dx) — W(x), A + — —(d/dx) — W(x), from which immediately follows that 
T is self-adjoint, whereas the SUSY algebra implies that the spectrum of T is 
non-negative. Being T a SUSY Hamiltonian, the Schrodinger equation l|40p can 
be explicitly solved. It turns out that the spectrum has a single discrete mode 
tpo for w = 0, and a continuum part ipk for u) > ujqCq [I]. 

The discrete mode corresponds to soliton translations, whereas the contin- 
uos modes are not reflected by the potential in Eq. (j40|) and describe soliton 
deformations. The set {^OjV'fc} is a complete, orthonormal basis and can be 
therefore used to expand a generic perturbation tp(x,t). 

Let us now consider our generalized SG model (j2"4")) . To simplify the calcula- 
tion we will first consider a vanishing coupling function F, later we will discuss 
the most general case. For F = 0, introducing the new dimensionless variables 
t' = uj c s t and x' = ccqX, with coq — K/b and the dimensionless parameter 
a = a/b, the field equations stemming from the Lagrangian (|24[) are 



l_dG 



c 2 \ 1 dG 

l + a^\ 6 vt > - (1 + a) 9 X , X , = - sin d + — — . (41) 
c t J u Q dV 

We want to investigate the classical stability of the 8 = <& soliton solution of 
the previous equations, which propagates at fixed speed v = c t - It is therefore 
appropriate to choose a frame in which this soliton is at rest. This can be easily 
done by performing in Eqs. (|4ip a boost of speed Ct, 



y = 1 (x' + %'), T = 7 (t' + ^x'), 7= , 1 (42) 

^/l-(f) 2 

In the (r, y) frame, Eqs. (|4"Tj) read 

<5>tt - $yy = -Sin$ + 2 

aa CJq a<i> 

/ „2 \ _ 1 

^rr - yy + a 1 + -§ TT + 2a-e T e v = -smd + — — . (43) 

Notice that the third and fourth term of the l.h.s. of the second equation are 
responsible for the breaking of Lorentz symmetry. Seen in the (r, y) frame, the 
solitonic solution of fixed speed Ct is static; it is given by 

$o (y ) = do {y) = arctan e v . (44) 
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The classical stability of the previous solution is investigated by considering 
small perturbations a, x around the static solution, 

Q(v,T)=0o(y)+a(y,T), 6(y,r) = 6 (y) + x(v,t). (45) 

At first order in the perturbation Eqs. give, 

<r TT - <7 yy = -(cos6»o)cr + B(a - \) 

Xrr - Xyy + a (l + Xrr + 2a—6 T (9 ) y = - sin 9 - B{a - x), (46) 

where B is a constant given by B = (l/w§)(<9 2 G/9$ 2 )(0). 

For generic values of the parameter a, the system (|46[) is very hard to solve. 
For a << 1 we can treat the terms proportional to a as a perturbation. At 
zeroth order in the perturbative expansion we can neglect them. In this case we 
can write Eqs. (|46|) in operator form, introducing the Schrodinger-like operator 



R 





«--£( 




( X J 







(47) 



The operator R can be diagonalized by the linear transformation ipi = a — \, 
^2 = cr + x 

R=(-& + ™ al ° , V ,48, 

where the potentials V\ and V% are 

Vx = (1 - 2sech 2 ?/), V 2 = (1 + 25 - 2sech 2 y). (49) 

Separating the variables in the perturbations, i.e setting Lpi(y,T) — ipi(y)e llAJ1T , 
¥>2(y,T) = Lp2{y)e. lul2T , equation (|47p becomes a time-independent Schrodinger- 
like equation. 

There are two particular values of the parameter B which correspond to 
SUSY models, B = and B — 3/2. In both cases the operator R is self-adjoint 
and non-negative, implying stability of the background solution. The B = case 
corresponds to the minimal, decoupled case discussed in Sect. El The diagonal 
elements R\ and R2 of the operator R of Eq. (|48|) corresponding respectively to 
the perturbations tpi,ip2, are equal and coincide with the previously discussed 
operator T of Eq. (140]) . The superpotential is also the same and is given by 
W% = W2 = tanhy. Ri and R2 have therefore the same spectrum, which 
coincides with that of the operator T. 

For B = 3/2, the operators R\ and R2, although not equal, correspond both 
to SUSY Hamiltonians. R\ is equal to the operator T of Eq. (|4T)]) whereas 
i?2 correspond to a different SUSY Hamiltonian (see e.g. Ref. [35j [36] ) . The 
superpotential is given by W2 = 2tanh(y). The spectrum of R2 has two discrete 
modes with eigenvalues UJ2 = and L02 — (3/4)wo c s- The continuum part of the 
spectrum is located at 0J2 > Aujqc s . 
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For generic B the operator Ri, being independent from B, corresponds al- 
ways to the SUSY Hamiltonian discussed above. On the other hand the operator 
i?2 does not seem to correspond to a SUSY Hamiltonian. From the general form 
of the potential V2, see Eq. . one can easily infer that we have a discrete 
spectrum for ui2 < (1 + 2B)luqc s , whereas we have a continuum spectrum for 
^2 > (1 + 2B)ojqc s . However, we do not have any argument to show that the 
eigenvalues of the discrete spectrum of R2 are non negative. In the following 
we will consider only the cases B = 0, but our considerations can be easily 
extended to the case B = 3/2. 

Let us now consider the first order in the perturbation expansion in the 
parameter a, i.e let us consider the system in the complete form P5|) . For 
B = the first equation of the system can be written as 

a TT — (T yy = —(1 — 2sech 2 (y)er. (50) 

Comparing with Eq. (|40j) . we immediately realize that this corresponds to the 
operator T = R\. It follows immediately that our fixed speed solitonic solution 
JH]) is stable for small perturbation of the field $. After some manipulations 
the equation for the perturbation x i n (|46[) (and B = 0) can be written as 

CXtt ~ Xyy + D(sechy)x T + (1 - 2sech 2 y)x = 0, (51) 

where C = l + a[l + (c s /ct) 2 ], D = a(c s /ct). The previous equation differs from 
the equation (|50p in two points. First, the term \ TT is multiplied by C . This is 
rather harmless, it just means that the eigenvalues of the operator T have to be 
rescaled by the positive factor C. This term cannot change the semi-positivity 
of the eigenvalues. Second, there is an y-dependent dissipative term. This term 
becomes negligible for y — > ±00. In fact in this limit Eq. (fSTj) gives 

CXtt ~ Xyy +X = 0, 

which has normal mode solutions x = Ae ltykyJrUT ^ with dispersion relation Cuj 2 — 
k 2 + l. 

The damping term is maximized near the minimum of the potential V(y), 
y = 0. Near y — the perturbation x nas a damping term \ oc e~' 3T , where 
(3 is some constant. It is interesting to notice in this context that the Lorentz 
symmetry breaking terms act on the perturbation as damping term and as 
rescaling of the eigenvalues of the operator T. 

The equation (|5ip can be solved by expanding x m terms of the eigenfunc- 
tions {^0,^} °f t ne operator T which represent a complete orthonormal set, 

/oo 
My)9k(T)dk. (52) 

Inserting this equation into (|5ip . using the orthonormal character relations, 
and the fact that {^0, V'fe} are eigenfunction of the operator T corresponding to 
eigenvalues and w 2 , one obtains after some manipulations the equations for h 
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and g/~] these turn out to be 



d 2 h n dh n 

+dJ°° dkdy[sech(y)^ k ^} + lu 2 9i = 0. (53) 

The equation for h is immediately solved to give 

h(r) = Pe- (3T + Q, (54) 

where P, Q are integration constants and (3 — 2ttD/C. From Eq. (JM]) it follows 
that after a transition time, only the pure translation mode ipo of the SG soliton 
will survive (when Q ^ 0); or that this mode will be completely damped to zero 
and only the static background SG solution will survive (for Q = 0). In both 
cases, this shows stability of our fixed speed soliton under perturbations of the 
discrete spectrum. 

The equation for the functions gi(r) in Eq. (|53[) is more involved. It is 
a integro-diffcrential equation that cannot be solved in closed form. However, 
from the general structure of the differential equation, it is evident that the 
solutions will always have an oscillatory part and a damping term ~ e~@ lT . 

Qualitatively the behavior of g will be therefore similar to the previously 
described one for h. After a transition time, the normal modes corresponding 
to the continuum spectrum of the operator T will be either completely damped 
or will survive as simple soliton deformations. This provides strong arguments 
for the stability of our solitonic solution for perturbations of the continuum 
spectrum. 

Until now we have considered explicitly stability when the coupling function 
F vanishes. Our discussion can be easily extended to the case when F ^ 0. 
Looking at the Lagrangian (j24|) one easily realizes that the terms introduced in 
this way are invariant under Lorentz transformations T(c s ). In the perturbation 
approach these terms have two effects. They change the value of the parameter 
C in Eq. (|5ip . and introduce a term proportional to sech 2 j/(cr — x) in both Eqs. 
(|46|) . The modification of C does not change the qualitative outcome of our 
previous discussion, whereas the modified Eqs. (|46|) can be analyzed using the 
same method we have previously used. 



9 Conclusions 

In this paper we have proposed a mechanism for fixing the velocity of relativistic 
solitons. The proposal has been first elaborated for a molecular chain model of 
double coupled pendulums. It has been then generalized to a full class of 2D 
field theories of the sine-Gordon type. 

From a phenomenological point of view, the model allows one to select the 
speed of a SG soliton just by tuning the elastic coupling constants and the 
kinematical parameters that support the transverse sound wave. Thus, the 
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selection speed mechanism could be at work whenever SG solitons appear with 
a given speed - and could be useful to formulate models whenever one needs 
to explain phenomena of this type, or however where the soliton speed plays a 
special role and should be tuned. Moreover, our selection mechanism could be 
also used to devise and realize non-linear media where SG solitons travel at a 
given fixed speed. 

Our results are also relevant from a fundamental, field-theoretical point of 
view. In fact we have shown that the speed selection mechanism is deeply re- 
lated both with the existence of some conditionally conserved quantity and to 
the breaking of the Lorentz symmetry of the usual SG models. The most strik- 
ing and highly non trivial result we have obtained in this context is that the 
characterizing features of relativistic SG solitons - that is, finiteness and local- 
ization of energy, topological classification and stability - may be still preserved 
even if the Lorentz symmetry is broken and a soliton of a given speed is selected. 

We have shown that the solitonic excitation with fixed speed of our model 
are classically stable and have topological conserved charges. In particular it 
turns out that the Lorentz symmetry breaking terms act on the perturbations 
around the fixed speed soliton as a damping term. 

This result is particularly intriguing in view of the analogy between rela- 
tivistic particles and SG solitons. It can be seen as a strong indication that 
relativistic nature of the particle is not a necessary condition for the analogy to 
hold. In particular one can speculate about the existence of SG solitons which 
behave as non-relativistic (e.g Newtonian) particles. 
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